Microscopic Theory of the Single Impurity Surface Kondo Resonance 
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We develop a microscopic theory of the single impurity Kondo effect on a metallic surface. We 
calculate the hybridization energies for the Anderson Hamiltonian of a magnetic impurity interacting 
with surface and bulk states and show that, contrary to the Kondo effect of an impurity in the bulk, 
the hybridization matrix elements are strongly dependent on the momentum around the Fermi 
surface. Furthermore, by calculating the tunneling conductance of a scanning tunneling microscope 
(STM) , we show that when the magnetic impurity is located at a surface the Kondo effect can occur 
with equal strength between bulk and surface states. We compare our results with recent experiments 
of Co impurities in Cu(lll) and Cu(lOO) surfaces and find good quantitative agreement. 
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I. INTRODUCTION 



When a magnetic impurity is located in the bulk of 
a metal it undergoes a non-trivial many-body scatter- 
ing with the electronic states at the Fermi energy, ep, 
called the Kondo efFecti. The bulk Kondo effect (to be 
contrasted with the surface Kondo effect studied in this 
paper) is one of the best studied problems in condensed 
matter physics and many different techniques from renor- 
malization group^ to Bethe ansata^ have been used over 
the years. The basic mechanism of the Kondo effect is 
the magnetic screening of the impurity (the formation of 
the Kondo singlet) at temperatures T below the Kondo 
temperature Tk- Above Tk the magnetic impurity be- 
haves paramagnetically but for T < Tk a resonance is 
formed close to the Fermi surface. The Kondo effect 
is very important in modern condensed matter theory 
and appears in many different areas of research, from 
U and Ce intermetallics (heavy-fermions)** to quantum 
dots^. There is a wide variety of phenomena that can 
be described within the universality of the Kondo prob- 
lem. The Kondo effect can be observed experimentally 
in measurements of the temperature dependence of re- 
sistivity (the so-called resistivity minimum)®, and also in 
thermodynamic measurements such as specific heat and 
magnetic susceptibility^ of dilute magnetic alloys, due 
to the enhancement of the density of states close to the 
Fermi energy (the Abrikosov-Suhl resonanca^). The en- 
hancement of the density of states is related to the change 
of the characteristic energy scales from e^^ to ksTK and 
therefore in the density of states from A^(0) (x 1/ep to 
7V*(0) cx l/{kBTK) > iV(0) since ep > keTK- 

The Kondo effect was also observed recently in STM 
studies of magnetic atoms on metallic surfacesSiSiiS (see 
Fig^. In a STM experiment electrons from a sharp tip 
tunnel into the material to be studied, creating a tun- 
neling current / due to the application of a potential 
V . Roughly speaking a STM experiment measures the 
local density of states at the Fermi energy via the dif- 
ferential tunneling conductance dl/dV—. When a STM 



tip is away from the magnetic impurity it measures the 
substrate density of states, iV(0); however, close to the 
impurity (that we call the adatom) electrons from the 
tip can tunnel directly to the impurity. The theory of 
STM is far from trivial because electrons from the tip 
not only tunnel to the adatom but also to the bulk and 
surface states, that is, there are various different channels 
of tunneling that lead to interference effects that have to 
be taken into consideration for the proper interpretation 
of the experimental data. 



STM Tip 
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Figure 1: Schematic picture of a STM measurement: zo is the 
distance of the adatom to the surface, R± is the distance of 
the STM tip from the surface, and is the distance, along 
the surface, between the STM tip and the adatom. 

The basic starting point for the study of the Kondo 
effect is the Anderson impurity Hamiltonian^^: 



J2 «fecLck<T 




(1) 



where Caa (clo-) is the annihilation (creation) operator for 
an electron on a localized atomic state (impurity state) 
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with angular momentum labelled by a, energy and spin 
a (cr =t, i), Ckcr (Ckjr) is the annihilation (creation) oper- 
ator for an electron on the conduction band with disper- 
sion Efc and momentum k, = c\„Caa is the number op- 
erator, U is the Coulomb energy for the double occupancy 
of the impurity state, and ^ka is the hybridization energy 
between the impurity and conduction states (the electron 
operators obey anti-commutation rules: {ca,c|^} = 5a,b)- 

One of the characteristics of the Anderson impurity 
model is the distinction between substrate and adatom 
wavefunctions. Although most theoretical works do not 
question the distinction made a priori between these 
quantum states, it turns out that this distinction is not 
completely natural. The reason is very simple: when an 
impurity atom is introduced in a metal, it hybridizes with 
the metallic states losing its identity. However, it leaves 
behind a phase shift 5k in the original metallic states. 
Thus, the impurity state cannot be really distinguished 
from the host states from the quantum mechanical point 
of view. Nevertheless, Andersoni^ makes the point that 
because the d-orbitals are a inner shell, the Coulomb en- 
ergy U for double occupancy of those orbitals is large and 
they must be distinguished from metallic states where the 
Coulomb energy is strongly suppressed by screening ef- 
fects. Thus, the distinction between these two types of 
states can only be clearly made when these states are 
orthogonal to each other so that the impurity does not 
cause a direct perturbation in the substrate spectrum. 
Nevertheless, even in metals where the electronic bands 
are generated out of d-orbitals (as in the case of Cu), 
the strong metallic bonding leads to a large s-wave char- 
acter of the bands and to very small overlap with the 
adatom d-orbitals (that is, these states are "naturally" 
orthogonaliSi) . Nevertheless, this orthogonality can only 
be distinguished a posteriori. In fact, we have found, 
by direct numerical computation, that this is the case 
in the surface Kondo effect. Finally, as pointed out by 
Andersoni^ (see also ref^AiL), the orthogonality of these 
states is not fundamental for the magnetic phenomenon 
which is essentially a local effect and all the subtleties as- 
sociated with orthogonality become encapsulated into the 
hybridization matrix elements ika which become a phe- 
nomenological parameter to be obtained indirectly from 
the experiment. However, in trying to understand the 
STM experiments, and especially the role played by the 
surface and bulk states, we cannot simply take these ma- 
trix elements as phenomenological parameters since we 
would not be able to separate the contributions coming 
from the bulk and the surface of the Cu substrate. 

Hence, in this paper we follow the scheme pro- 
posed by Tsvelik and Wiegmannii and work with non- 
orthogonalized states. In this case the hybridization en- 
ergies are given by: 

<ka = (efe + ea) J d\ V':(r)Vk(r) (2) 

where V'a(r) is the atomic state of the impurity which 
is localized on the impurity site (that is, it decays ex- 



ponentially away from the impurity) and V'k(r) is the 
Bloch wavefunction of the crystal. The hybridization en- 
ergy depends not only on the amplitude but also on the 
direction of k. Nevertheless, it is usually assumed that 
the hybridization is constant in momentum space, that 
is, ik,a = t- Although this assumption maybe warranted 
for spherical gapless Fermi surfaces (since the momen- 
tum dependence "averages out" and most of the interest- 
ing physics occurs at the Fermi surface) it is not valid 
for metals where the chemical potential crosses gaps in 
the electronic spectrum. The Kondo Hamiltonian can 
be obtained from in the limit of \t\ <C Ca, C/ via the 
Schrieffer-Wolf transformation^^ that maps the Hamilto- 
nian onto: 

Hk - efecLck. + Jk^- s(0) (3) 

where S is the impurity spin, s''(0) = 
V2Z]k,k',a,/3Cka'''a,/3'^k'/3 is the electron spin at the 
impurity site (cr" with a = x,y, z are Pauli matrices) 
and 

Jk « If + (4) 

\\ea\ \ea + U\J 

is the exchange interaction between local spins and con- 
duction electrons (the energy is measured relative to the 
Fermi energy). In going from to Q we have disre- 
garded a series of physical processes that cannot be de- 
scribed as spin exchange interactions, foremost being the 
variable occupancy of the impurity site. 

Even simple metals like Cu have energy gaps in their 
electronic spectrum due to the periodic potential gener- 
ated by the Cu ionsiSiii. For instance, it is well-known 
that for Cu along the (111) direction there is a gap of 
5.1 eV at the Fermi energy while there are no gaps along 
the (100) direction (see Figl^J. On the one hand, tka 
does not exist for bulk states when k points along the 
(111) direction. However, because of the bulk gap, one 
has a surface state in the Cu(lll) surfacei^ (a state that 
decays exponentially in the (111) direction but it is ex- 
tended in the directions perpendicular to (HI)) that can 
hybridize with the impurity located on the Cu(lll) sur- 
face. On the other hand, there are no surface states in the 
(100) direction to hybridize with an impurity located on 
the Cu(lOO) surface, although bulk states hybridize with 
it since their wavefunctions decay exponentially outside 
of the crystal. Previous studiesi^ on the surface Kondo 
effect either assume a momentum-independent hybridiza- 
tion tka or sum average over a gapless Fermi surfac&ii. 
By ignoring momentum dependence or the gap structure 
of the Fermi surface, the comparison between Cu(lll) 
and Cu(lOO) surface Kondo effects of noble metals is not 
possible. In fact, we show that recent STM experiments 
of atoms located in the Cu(lll) and Cu(lOO) surfaces 
have been interpreted incorrectly because of these differ- 
ences. In particular, we show that the surface Kondo 
effect between a magnetic atom located on the Cu(lll) 
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surface and the Cu(lll) surface state is quantitatively 
similar to the surface Kondo effect of a magnetic atom 
on a Cu(lOO) surface with bulk states. The main reason 
for this similarity comes from the fact that the bulk state 
that dominates the Kondo scattering is the one with the 
largest penetration in the work-function barrier, that is, 
the state with momentum perpendicular to the Cu(lOO) 
surface. The state with momentum perpendicular to the 
Cu(lOO) surface has essentially the same decay rate out 
of the crystal as the surface state in the Cu(lll) surface 
and therefore similar hybridization with the magnetic im- 
purity. 



^ K 



X r K 



Figure 2: Planar cut of the bulk Cu Fermi surface in the first 
Brillouin zone: F-L is the (111) direction (notice the absence 
of states due to the bulk gap); F-X is the (100) direction. 

In the presence of the tip new terms have to be added 
to the Hamiltonian (1) that describe the tunneling pro- 
cesses of tip-to-adatom and tip-to-substrateii. The tun- 
neling processes of tip-to-adatom is given by: 

Jiat = ^ tap (cLcpa + H.c.) . (5) 

(T,a 

where c^a (c],^.) annihilates (creates) and electron at the 
tip and tap is the hybridization energy between tip and 
adatom. The tip-to-substrate hybridization is given by: 

i^i = 5I^kp(4.Cpa + H.c.), (6) 

ka 

where tkp is the hybridization energy between substrate 
(bulk or surface state) and adatom. The tip Hamiltonian 
is simply 

a St = XI ^V^loCpa (7) 
a 

where ep is the energy of the electron states on the tip. 
The total Hamiltonian of the tip-substrate-adsorbate sys- 
tem is H — Hg + Ht + Hat + Hst ■ As in the case of (Q) 
the hybridization energies are given by: 

W = (ea +e/3) J dVV';(r)V'/3(r) (8) 

with a, P = k, a, p labeling the wavefunctions, V'a(r), for 
substrate, adatom and tip, respectively. Notice that the 



complexity in this problem comes from the interference 
between the different channels of tunneling. 

In the present work, we will consider the case of Co on 
Cu surfaces while our theory is applicable to any mag- 
netic adatom on noble metal surfaces. Our main objec- 
tives are: (i) the calculation of the hybridization ma- 
trix elements between the different wavefunctions in the 
problem (bulk, surface, impurity, and STM tip) with the 
smallest number of adjustable parameters, and (ii) the 
description of the surface Kondo effect via STM mea- 
surements. While the spectrum of bulk and surface states 
is readily available in the literature, very little has been 
published on the actual form of the states. Instead of em- 
barking on a complicated calculation of wavefunctions via 
heavy numerical techniques, we opted for a simpler ap- 
proach that can provide quantitative results that can be 
directly compared to the experiment as well as intuitive 
understanding of the problem. As we are going to show, 
our simplified models have limitations in explaining some 
of the available STM data and more detailed work is re- 
quired. We hope our work stimulates other works that 
make use of more sophisticated methods to calculate the 
matrix elements because, as we are going to show, they 
are fundamental for the understanding of magnetism in 
metallic surfaces. We model the periodic potential inside 
of the crystal in the nearly-free electron approximation 
and obtain the bulk gaps as observed in metallic Cu. 
The surface states are obtained after the modeling of the 
surface potential in terms of the image charge in the vac- 
uum. The energies of surface states and resonances in 
each surface are in very good agreement with photoemis- 
sion results in Cu(lll) and Cu(lOO) surfaces. Following 
the literature in STM theory, the impurity is modeled via 
standard Hydrogen-like wave functions. 

The paper is organized as follows: In Sec.m we model 
the ionic potentials that give rise to the wavefunctions 
of the metallic host, adatom, and the STM tip. Using 
the wavefunctions and their spectrum we discuss the Cu 
band structure, both the bulk and Cu(lll) surface state. 
In Sec. mil we present the results for the hybridization 
matrix elements that determine the hybridization ener- 
gies as a function of the distance from the surface and 
adatom and also as a function of the momentum. The 
surface Kondo effect is discussed in Sec. lIVl in the context 
of the Schrieffer-Wolf transformation and the slave-boson 
mean-field theory. Using the calculated hybridization en- 
ergies the differential conductance is calculated in Sec.lVl 
and compared with the STM data. Sec. I VII contains our 
conclusions. 



II. WAVEFUNCTIONS AND SPECTRUM 

In this section, we obtain wavefunctions needed to cal- 
culate the tunneling matrix elements This is the 
first step towards the STM differential conductance. As 
discussed in the introduction we are going to use simple 
models for the various potentials involved in the problem 
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so that we can make analytical progress in calculating 
the hybridization energies. 

Metallic Cu is a very good metal with a nearly spherical 
Fermi surface except for the gaps in the (111) direction. 
These gaps have their origin in the periodic potential gen- 
erated by the ions (see FiglSJ. Although the Cu states can 
be obtained from first principle calculations-'^^ we adopt 
the nearly free electron approximation and model the 
bulk potential with a simple periodic function. Further- 
more, in order to study the surface states we consider a 
semi-infinite crystal where the bulk is located at z < 
(see Fig. 121. We write the potential in the substrate as: 



V{r) 



'2Vi cos 2^ 
-2Vi 



Vo 



4(z- 



, for z < 0, 

, for < z < zi, 

, for z > zi, 



(9) 



where Vi is the bulk potential, a is the lattice spacing 
in the direction of the surface, Vq is the work function 
measured from the bottom of the conduction band, Zim 
is the image plane, chosen in such a way that the image 
potential joins the constant potential —2Vi continuously, 
and is related to zi as 



4(^0 + 2^1) 



(10) 



In the potential Eq. ((HJ, zi is the only free parameter. 
The eigenfunctions and eigenenergies for the Schrodinger 
equation associated with the potential ® can be calcu- 
lated exactly and the details are given in Appendix A. 



provide a very good description of the band gap in the 
(111) direction and also of the ionization energy of metal- 
lic Cu. The Cu(lll) surface has a band of surface states 
as measured by photoemission, starting 390 meV below 
the Fermi surface- By solving the SchrSdinger prob- 
lem for the potential © we find the surface band (see 

Appendix A) and fix the value of zi — 1.66 A in order 
to fit the photoemission results. Using this value for zi, 
we also obtain the next surface state at 788 meV below 
the vacuum level. This result agrees well with photoe- 
mission experiments that find the second surface state at 
830 meV below the vacuum levelSi. Thus, we were able 
to get a consistent description of two independent exper- 
iments with a single fitting parameter. In Fig. ^we show 
a plot of the surface state wavefunction as a function of 
the distance away from the surface (the wavefunction is 
chosen to be a real function). Notice the slow decay of 
this state inside the crystal in the direction perpendicular 
to the surface. 
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Figure 3: Plot of the image potential energy near a crystal 
surface. The solid dots are crystal atoms. 

Let us consider first the case of the (111) direction. 
In order to fit the experimental band structureii along 
the (111) direction in the nearly- free-electron potential 
we find that Vi = 2.55 eV, Vq = 13.55 eV. These values 



Figure 4: The Cu(lll) surface state with energy 0.39 eV be- 
low the Fermi surface. The vertical axis represents its nor- 
malized wavefunction multiplied by the box size L. The hor- 

o 

izontal axis represents z in A . 

In the 3-dimensional fc-space of Cu, a gaped "cavity" is 
centered around the (111) direction near the Fermi level 
and has the bulk band surrounding this gaped region^-. 
By decomposing the momentum of a bulk state k into 
k_L (perpendicular to the surface) and ky (parallel to 
the surface), one expects that the larger k±, the larger 
wavefunction tail outside the crystal, which makes a sig- 
nificant difference on their contributions to the surface 
Kondo effect. As a result, it is important to carefully ex- 
clude the (111) gaped region from the Fermi sphere when 
calculating the surface Kondo effect. 

It is known that a nearly-free-electron band has a 
parabolic dispersion not too close to the band edge, where 
"band edge" refers to the interface between the bulk 
band and the gap in the 3D fc-space. When calculat- 
ing the STM tunneling conductance, one would have to 
include the contribution close to the band edge as well. 
However, the band-edge contribution is typically much 
smaller than the one coming from the parabolic band, so 
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we can take merely the contribution from the parabolic 
band. For the Cu(lll) band, we terminate the parabolic 
band 0.7 eV below the band edge. As long as we work in 
a parabolic band, the bulk-state wavefunction inside the 
crystal takes a plane-wave form. 

We now turn to the case of the (100) direction. We 
find that we get good agreement with band structure 
calculationsiZi by choosing Vi = 3.05 eV, and Vq = 13.45 
eV. In contrast to the (111) direction, the (100) direc- 
tion has no true surface state near the Fermi surface. 
However, there is a surface resonance at 1.15 eV above 
the Fermi level22i that can be well described by choosing 

o 

zi = 1.51 A . When calculating the tunneling conduc- 
tance, we can safely sum over the entire Fermi sphere and 
adopt a parabolic dispersion, i.e. a plane-wave wavefunc- 
tion inside the crystal. 

It is the d-orbital of the adatom that actually partici- 
pates in the Kondo resonance. When one of these atoms 
is placed on a metal surface, its outmost s-wave electrons 
either get transfered to the metal conduction band or to 
its own d-orbital. For example, the Co adatom tends to 
form [ArjBd^ on a Cu surface^^ and therefore has an ef- 
fective spin 1/2. As a result, an adatom d-wavefunction 
can be modeled as a Hydrogen atom with an effective 
charge ZeS given by: 

V'a(r) = i?„2(^off, r)Y2,n{e, 0) , (11) 

where Y2m{d^(j)) are Spherical Harmonics and i?„2 are 
Hydrogen radial wavefunctions. The effective charge Z^s 
can be determined by comparing rmax, defined by, 

dm'' 



dr 



= 0, 



(12) 



with the experimentally observed d-orbital radius for 
metallic Co, r^ax = 1-25 A. In our calculation we ori- 
ent the axis of quantization for the atomic problem along 
the direction perpendicular to the surface. 

The STM tip wavefunction is modeled following Tersoff 
and Hamann^** who proposed a wavefunction of the form: 

^K.R exp(-Kr) 

(13) 



1, 



for r > R, 
for r < R, 



where R is the radius of curvature of the tip about its cen- 
ter, the decay constant k = •\/2mt (0 — e^) controls the 
wavefunction tail, and is the tip work function. Notice 
that the adatom and tip wavefunctions are normalizable, 
that is, / d'^r|i/'a(r)p = / d^r|V'p(r)p = 1, when inte- 
grated over the entire space. 

The set of three different wavefunctions (substrate, 
adatom, and tip) are the main elements required for the 
calculation of the hybridization energies in the Anderson 
impurity Hamiltonian. 



III. HYBRIDIZATION ENERGIES 

Using the results of Eq. (|A1(I , Hll|) and (|13|) we can eval- 
uate the overlap integrals in © that are important for the 



hybridization energies. Let us first remark that the nor- 
malization of each of these wavefunctions is different; for 
instance, the bulk state is normalized in a box of size L^, 
where L is one of the dimensions of the box. The surface 
state, however, since it decays exponentially away from 
the surface but is extended on the surface, is normalized 
in an area of size L^. Finally, the adatom and tip states 
always decay exponentially and therefore do not require 
any box normalization. Thus, it is convenient to redefine 
the wavefunctions so that: 



1 



L3/2 



03,k(r) 



V'S,k(r) = -02,k(r) 



(14) 



where the subscript B refers to bulk and S to surface. 
Note that in this case the hybridization energies between 
substrate (bulk and surface) and adatom in (jSJ can be 
written as: 



tk,a.a = (ea+e(k))-^l/i'^)(k) 



(15) 



where d — 3 for bulk and d — 2 for surface. The matrix 
elements are, therefore. 



l/i'^)(k) - / dV0rf,k(r)Va(r). 



(16) 



Notice that Va'^^ (k) is not only a function of the momen- 
tum and the distance zq between adatom and surface, 
but also depends on which surface - here labeled by the 
subscript a - the adatom is located. 

In Fig. Owe show the result for Vf^iii-) , the hybridization 
matrix element between m = adatom state and surface 
state in the Cu(lll) surface, as obtained by numerical 
integration of (|16|l using (|A1I) and Hill) , as a function of zo 
(in units of the distance between the planes in the (111) 
direction) by assuming k to be at the Fermi surface of 
the surface band [since the surface states have essentially 

(2) 

a spherical Fermi surface there is no dependence of V^il 
with the direction on the Cu(lll) surface]. The first 
striking result is the strong oscillation of the value of the 
matrix element with zq. This oscillation is the result of 
the interference between the d-state of the adatom and 
the surface state. Naturally, the largest overlap occurs 
at the surface (where the surface state is maximum, see 
Fig. 2)) but we notice that there is substantial overlap 
between the surface state and the adatom at one lattice 
spacing from the surface. In principle, we cannot tell 
what is the actual orientation of the adatom relative to 
the surface but we have checked that the m = state 
is the one with biggest overlap as compared with higher 
angular momentum states, thus, instead of a sum over the 
angular momentum states in Eq.JQ) we could have kept 
only the one with m — 0. Furthermore, first principle 
calculations for adatoms in metallic surfaces indicate that 
the atomic orbitals are oriented in such a way to generate 
the largest overlaps^. 
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Figure 5: V/f^ dashed-dotted line, V/fi'(r-A'): thick line, 
and Vljl(r-X): thin line, as a function of . 



In Fig. we also plot the value of V'^^i) , the hybridiza- 
tion matrix element between the bulk states and adatom 
in the Cu(lll) surface as a function of the distance from 
the surface for two different bulk states, all of the them 
at the Cu Fermi surface, with momenta along the F-X, 
and F-K directions (see Figl^J- Once again we notice the 
oscillations due to the interference between those states 
and the adatom states but we also observe that the over- 
lap now is a factor of 3 smaller because these states do 
not have a large amplitude at the surface. Another inter- 
esting result of these calculations is the fact that as the 
momenta moves towards the (111) direction the overlap 
grows. This is due to the fact that away from the (111) di- 
rection the Fermi surface of Cu is essentially spherical and 
therefore the states at the Fermi surface have essentially 
the same amplitude, however, the states with large mo- 
mentum along the direction of the Cu(lll) surface have 
bigger penetration on the potential barrier generated by 
the work function and therefore larger hybridization with 
the adatom. This effect will be very important later for 
the interpretation of the experimental data. 

In Fig. we show the value of V'^joo) overlap 
between the state of an adatom on the Cu(lOO) surface 
with bulk states with momenta along the F-X (this direc- 
tion is perpendicular to the Cu(lOO) surface), and F-K 
(since there are no surface states for the Cu(lOO) surface 
we do not need to calculate y^^^ for this surface). Be- 
sides the oscillations found in the previous calculations, 
we find that when the momentum is oriented perpen- 
dicular to the surface its penetration is maximum and 
therefore has largest overlap. In fact, we find that the 
overlap between the adatom state on a Cu(lll) surface 
with the surface state is essentially the same value as the 
overlap between the state of the adatom on the Cu(lOO) 
surface with a bulk state with momentum oriented in 
the (100) direction. One can easily understand this ef- 
fect from the nearly free electron picture: in the absence 
of the periodic potential of the crystal, the Fermi surface 
states in the (111) and (100) direction would be essen- 



tially identical and their decay in the vacuum of their 
respective surfaces would be the same. When the peri- 
odic potential of the crystal is added a gap opens in the 
(111) direction but not on the (100) direction. However, 
the bulk state that disappeared at the Fermi surface in 
the (111) direction now becomes the surface state in this 
same surface. Therefore, its amplitude at the surface 
and close to it is essentially the same as the state in the 
(100) direction (this is can be readily understood due to 
the long decay rate of the surface state shown in Fig. 
and the overlap is essentially the same. Thus, one has 
to be careful when discussing what happens to adatoms 
in different surfaces because there is very little difference 
between surface states along gaped directions and bulk 
states along non-gaped directions in the Fermi surface of 
noble metals. 



2 3 
> 




-4 -2 2 4 



Figure 6: Vlo^{r-X): thin line, ^^^(r - K): thick line, as a 
function of zo- 



IV. THE SURFACE KONDO EFFECT 

There are two main limits of interest in the Kondo 
problem. The first one is reached when |ika| ^ |ea|,C^ 
and we can treat the hybridization energy as the smallest 
scale in the problem. This is the so-called Kondo regime 
where the average occupation of the impurity state is 
essentially constant, that is, (ua) « 1. In order to study 
this limit it is convenient to split the electron states into 
bulk (B) and surface (S) so that we write the Anderson 
impurity Hamiltonian as: 

H = Xl^lWcI,k^Cl^k<T +X1^2(q)4_q^C2,q<. 
k,(T qo- 

-I- eg dldg- + U ng,! rig^x 

a 

k,o- 

+ 5I*'S.q-(4,qarf- + ^-c-) (17) 

qa 
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where ci^ko- (c2,qtT) describes the bulk (surface) electron 
annihilation operator, respectively. We have a two-band 
problem where both bulk and surface states hybridize 
with the impurity at the surface. Notice that this is a 
single channel problem since only one effective channel 
couples to the impurity. 

In the Kondo limit we can use a Schrieffer-Wolf 
transformation^^ and rewrite the Anderson impurity 
Hamiltonian as a exchange problem: 

H ^ en(k)c,\ k,TC",k^ 

k,(T,ra=l,2 

+ E E S-aa,;3clk,„^n,m(k,k')c„,a^) 
k,k' a,l3=-\.,l n,m=l,2 

where S is the localized impurity spin and 



<Ai,m(k, k ) — 



1 



^(d„)(k)i/(d™)(k') 



X 



(19) 



is the exchange energy matrix (we have defined di — 3 
and ^2 = 2 as the dimensionalities of bulk and surface 
states, respectively), and 



1 



1 



U 



(20) 



Notice that, unlike the single band Kondo effect elec- 
trons from one band can scatter from the impurity be- 
tween two different bands with spin flip (all energies are 
measured relative to the Fermi energy). If we ignore the 
off-diagonal matrix elements of p9|l the problem trivially 
reduces to two independent Kondo effects with character- 
istic Kondo temperatures 



-1/9- 



(21) 



where Z)„ is the bandwidth cut-off (Di 
0.4 eVi'') and 



.6 eV, D2 



gn 



1 KMlV^'-HkrW 



X 



(22) 



is the effective coupling constant, Nn{0) is the Fermi sur- 
face density of states (for spherical Fermi surface we have 
iV„(0) = (L/7r)'^"27rm*A:^""Vfi^ where m* is the effective 
mass). Because of the difference in the Kondo tempera- 
tures given in H21|l one of the scattering channels domi- 
nates the screening process and the singlet state can be 
either formed with the surface state when Tk,2 > Tk,i 
or with the bulk states if Tk,i > Tk.2- If we apply this 
argument to the results of the overlap integrals of the pre- 
vious section we see that for the adatom in the Cu(lll) 
surface the Kondo screening is dominated by the sur- 
face state since the coupling constant (E^ is bigger for 
this state than the bulk states. On the other hand, if 
the adatom is sitting at the Cu(lOO) surface, the Kondo 
effect is dominated by the bulk state with momentum 
perpendicular to the Cu(lOO) surface. 



The above discussion is valid as far as the occupation 
of the impurity level is close to 1. However, in the prob- 
lem at hand the adatom occupation can change substan- 
tially with the hybridization since the matrix element is 
strongly dependent of the distance to the surface and 
the type of states involved, that is, tka can be of the or- 
der of ea (but still much smaller than U). In this case, 
the mapping into the Kondo problem fails to properly 
describe the actual occupation on the adatom. Instead, 
we should start from the Anderson impurity Hamiltonian 
([T|) and compute the Kondo temperature directly from it. 
Here we will do so by adopting the large N calculation of 
Newns and Read2& that allows for average fluctuations in 
the valence of the adatom in the limit of t/ ^ 00 (in fact, 
estimates based on LSDAiS give a value of C/ « 2.84 eV 
and ta ~ —0.84 eV). The idea of the large N calculation 
is to generalize to include more degenerate electron 
states (besides the 2 states associated with the spin 1/2) 
and take the limit where the degeneracy diverges. Math- 
ematically this is accomplished by generalizing the spin 
index in the electron operator to a new index m that 
varies from 1 to N . It is convenient to write Anderson 
impurity Hamiltonian (^) in the limit of [/ ^00 as: 

m k,m 

P + H.c.) (23) 

k,m 

where V is the projector operator that projects out any 
states of the adatom with double occupancy. In the 
adatom we can, in principle, have three different states 
with zero, \dP), single, and double, \(P), occupancy. 
We now expand the Hilbert space by introducing a "phan- 
tom" boson state such that when the state of the adatom 
is unoccupied by an electron it is occupied by a boson, 
say and when it is unoccupied by a boson it is 

occupied by an electron, \}pd}). Observe that this con- 
struction constrains the number of bosons and electrons 
to be 1 for all states IV-') in the Hilbert state: 



{na + b^b)\^l:) = 1| 



(24) 



The utility of this representation becomes clear when we 
consider the action of the operator c\^h on the state 
\b^dP) producing the state \lPd}). Further application of 
c\jyjD to \lPd}) annihilates the state while application of 
Cambi^ to \b^d}) returns the system to the state \b^<P). 
Thus, the repetitive application of cjj„& or Camb'' never 
produces the doubly occupied state. In other words, we 
can replace V in (|23|l by a boson operator if we ensure 
that the constraint (|24|l is obeyed. This constraint can be 
seen as an energy cost, that is, states of the Hamiltonian 
that violate the constraint (|24|) are highly excited states. 
If we are interested in the low energy physics then we 
can convert (|24|l into a new energy scale and add it to 
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the Hamiltonian as: 

m k,m 
+ X! ^^-^ ("^km^am^ + H.C. 

+ 7K + 6^6-1) 



where 



(25) 



where 7 is a Lagrange multipher. Notice that in the hmit 
of 7 — > CX3 the constraint (|24|l is automatically obeyed 
since the energy is minimized when (I24II is fulfilled. As 
is, Hamiltonian H25() is still very complicated and at this 
point we have to use the large N approximation that 
consists in treating H25|l in a mean-field approximation. 
In the mean-field limit we consider the average value of 
the boson fields instead of their actual value, that is, we 
replace the operators b and by 



(6) = (6t) 



(26) 



Notice that by treating the bosons in the mean-field ap- 
proximation we obtain a quadratic Hamiltonian in H25() 
that can be diagonalized by standard techniques and the 
total energy of the system can be minimized with respect 
to z and 7. The physical meaning of these parameters be- 
comes clear when we calculate physical observables. For 
instance, the average occupation of the impurity is given 
by: 



= I - z 



(27) 



that shows that z has to do with the change in the occu- 
pancy of the adatom state due to the hybridization with 
the substrate. On the other hand we can also show that 7 
is the renormalization of the energy of the impurity state, 
that is, the renormalized impurity state energy is given 
by: 

— 

+ 7. The mean-field parameters Ca and z 
are calculated by solving a set of coupled transcendental 
equations^^: 



2 {zFakaiep] 

1 — z — — arctan — — 



:-F'afca(eF) In 



V 



{4 + [zFaU^F)f} 

D 



1/2 > 



where -F'Qfe/3(w) are matrices defined, in general, as: 



w — e^ + iri 



(29) 



and D is bandwidth of the conduction band (notice that 
a,/? = a,p in reference to the adatom and tip, respec- 
tively). The Kondo temperature can be directly ex- 
tracted from this formalism aaS^ 



ksTx = zAo = L'cxp • 



(30) 



Ao = liaFakaiuj). 



(31) 



Therefore, using the large N account we can deal with 
the case where the hybridization energies are of the order 
of the adatom energies. We use these results in studying 
the STM differential conductance in the next section. 




0.25 0.5 0.75 



1.25 1.5 



Figure 7: Kondo temperatures as functions of the Co distance 
from the surface, zq, on the Cu(lll) surface. 



We are not aware of any measurement of Co-to-Cu dis- 
tance being reported, so this distance zq will be treated 
as a parameter of the order of the Cu lattice space. 
Fig. [7| and |H1 shows our calculated Kondo temperature 
as a function of zq- An interesting feature is that the 
Kondo temperature decays fast with zq, which can be 
understood because of the exponential overlap between 
the substrate and Co wavefunctions. When calculating 
the Kondo temperature for Co/Cu(lll) assuming that 
only the bulk state exists, we find Tk < lO^^K for the 
zq's as Fig.[7| which clearly means that the surface state 
dominates the Kondo effect. The domination is not sur- 
prising since a significant portion of large-fc_L bulk states 
are gaped around the (111) direction. Comparing the Cu 
Fermi surface with a perfectly spherical one, one finds 
that the Cu surface state contributes to the Kondo effect 
in a way similar to the bulk band of the spherical Fermi 
surface around the (111) direction. Similarly, the bulk 
band with respect to the Cu(IOO) surface can be consid- 
ered as a perfectly spherical Fermi surface for calculating 
the Kondo effect. As a result, the Kondo effect of the 
Co/Cu(lll), dominated by the surface state, should have 
the same order of magnitude as Co/Cu(100). The posi- 
tions that provide the experimentally measured Kondo 
temperatures2i (TV ^m) — 54K and T/f,(ioo) = 88K) are 

o 

0.41 A and from the (111) and (100) surfaces, respec- 
tively. Hence the STM conductance will be calculated by 
fixing Co atom at these optimized positions. 

We should point out that from an experimental point 
of view we would expect the Co adatom to be located 
around one lattice spacing away from the surface. Thus, 
our obtained values of their position seem to be quite 
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small (especially in the case of the (100) surface where 
the adatom has to be located at the surface) and therefore 
unrealistic. Nevertheless, one has to have in mind that 
the there are many assumptions made in the calculation 
that can affect significantly the value of the Kondo tem- 
perature, namely, the form of the image potential close 
to the surface, the large N approximation, and the limit 
of U ^ oo. Furthermore, it is clear from H3U|I . that the 
Kondo temperature is exponentially sensitive on the de- 
tails of the hybridization and energy scales and there- 
fore is significantly affected by the various assumptions. 
Thus, the plots in Fig[7|and Fig|Slhave to taken as the 
trends for the variation of the Kondo temperature with 
the distance from the surface and not as realistic values. 



80 - 



where the differential conductance with a subscript "0" 
refers to the background signal, A is the amplitude of 
the STM conductance. 



A cx {ImFpkaf (1 + g') , 
q is the Fano lineshape parameter, 

_ tpa + Rc Fpka{i^) 

and ^ is the dimensionless bias voltage 

eV + la 



(35) 



(36) 



(37) 



Expression has been frequently used for the fitting 
of STM data. 
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Figure 8: Kondo temperature as functions of the Co distance, 
zo, on the Cu(lOO) surface. 




V. STM DIFFERENTIAL CONDUCTANCE 

We assume that the STM tip and the adsorbate- 
substrate complex are each in local equilibrium, and use 
the standard Green's function formalism for calculation 
of the differential conductancaii: 

^ cx -Im {Fpkp{L0) + G'q(w) [tpa + Fpka{'^)f\ 

(32) 

where F^tpi'^) is defined in H29|) and Ga is formally the 
adsorbate retarded Green's function associated with the 
total Hamiltonian H. However, because the effect on the 
adsorbate from the tip is much smaller than that from 
the substrate, it is always safe to regard Ga as the one 
without tip. In this work, Ga (without tip) is calculated 
by the mean-field slave-boson technique2& as: 



Ga{uo) 



1 



(33) 



where the symbols here are given in the previous section. 

By substituting lO into the STM differential 

conductance can be written as: 



dl f dl , , 



q 



1 + 2q( 



(34) 



Figure 9: Differential conductance (in arbitrary units) at 
7?ll = 0, as a function of the bias V in meV for an adatom on 
the Cu(lll) (thick line) and Cu(lOO) (thin line) surfaces, as 
obtained from the theoretical calculations. 

The shape of the differential conductance versus the 
bias voltage is determined by the Fano line shape param- 
eter q. In FigEl we show the differential conductances 
obtained from our microscopic theory for the Cu(lll) 
and Cu(lOO) surfaces. From Eq. H36|l we see that q is re- 
lated to three microscopic quantities tpa, Im Fpkaiijj), and 
R-eFpfca(^)- It is found that tpa <C ReFpka{^) and there- 
fore we can set tpa = in (|36|l . Obtaining the first two 
is straightforward since they only depend on the physics 
at the Fermi energy. From H29(l . we have: 



luiFak/s (w) 



Try^JaktkfjSjuj 

k 



(38) 



Notice that the amplitude of the STM signal, (|35|l . de- 
pends essentially on this quantity and can be calculated 
reliably. In Fig^|we plot A{R\^) for R± = 6 Aand com- 
pare it with the experimental data^i. As one can see, 
the agreement between the theory and experiment and 
theory is quite good. 

The other quantity of interest, namely, Re-Fpfco(w), re- 
quires an integration over the entire Cu band. For the 
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Figure 10: Amplitude of the differential conductance (in arbi- 
trary units) for Co on Cu(lll) as a function of -Ry as compared 
with the experimental data—. 

Cu(lll) surface, the surface band has a free-electron dis- 
persion starting 390 meV below the Fermi level up to 
400 meV above it, where the upper limit is the inter- 
cept of the surface band to the bulk bands. We calculate 
Ke Fpkai^) for this surface band and determine the tun- 
neling conductances, H34|l . for different values of as 
shown in shown in Fig. ^2 
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Figure 11: Differential conductances for the Co on Cu(lll) 
for different values of iiy . From bottom to the top: iiy — 
0,2,4, 6,8, 10 A 

We notice that for short distances away from the 
adatom, the theory can describe very well the lineshape of 
the tunneling conductance. Nevertheless, as the distance 
of the STM tip to the adatom is increased, the lineshape 
becomes more symmetric and a positive Lorentzian is 
developed (see Fig lll() . The origin of this result can be 
easily understood: notice that in FigEHthe amplitude of 
the STM conductance becomes very small for large val- 
ues of This smallness can be tracked down to the 
small value of ImFpka in (|35l) indicating that the inter- 
ference between tip and adatom is fading away. However, 
one sees from ()36|l that this indicates that the Fano pa- 
rameter q is becoming very large and positive given rise 



to a positive Lorentzian. This feature is general of any 
theory of the STM conductance that starts with nearly 
free electron wavefunctions and cannot be avoided since 
it has its origins on the oscillating nature of the surface 
states density of states^^. This result is at odds with the 
experiment when the distance between tip and adatom 
becomes large . Instead of a positive Lorentzian, it is 
observed a negative Lorentzian ■^^ . 

In Fig^lwe plot the calculated values of q against the 
values used by the experimentalists in order to fit the 
data—. The discrepancy between theory and experiment 
is clear. More work will be required in order to check this 
issue. 
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R|| 

Figure 12: Plot of the Fano parameter q{R\\) for Co on 
Cu(lll) as compared with the values obtained by the ex- 
perimental fit—. 

In order to calculate Re Fpka (w) for an adatom located 
at the Cu(lOO) surface we need information of the whole 
Cu bulk band. The Cu bulk bands start at approxi- 
mately 10 eV below the Fermi level and disperse up to 
the vacuum level that is located at approximately at 4.8 
eV above the Fermi energy. At the Fermi energy the 
Cu band has dominant s-wave symmetry which gives the 
quasi-spherical shape of the Fermi sea. However, there 
are also Cu d-bands located from 2 to 5.5 eV below the 
Fermi level with very large density of statesA^. Using the 
nearly free electron approximation (s-wave) with a cut-off 
at the vacuum level we have calculated ReFpka{uj) and 
have found that it yields a value of q{R\\) that changes 
sign from positive to negative as R\\ increases, which is 
not the case observed in the experiment. In principle 
the quantity Re Fpkai^) must be calculated including the 
d-bands. Using a simple parabolic dispersion for the d- 
bands we find their contribution to be less than 0.1% of 
the total value of Rei^pfea(w) and therefore these bands 
are not able to fix the problem of the sign of q{R\\)- Thus, 
we have problems with the use of bulk states in our cal- 
culation. We believe that band structure effects are im- 
portant for the exact calculation of the surface Kondo 
effect due to bulk states. As far as we know there are no 
such detailed calculations. 

In order to remedy this problem we introduce a cut- 
off in the problem such that by integrating bulk states 
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Figure 13: Amplitude of the differential conductance (in arbi- 
trary units) for Co on Cu(lOO) as a function of -Ry as compared 
with the experimental data—. 

starting from the bottom of the nearly free electron band 
we stop the integration when the value of Re Fpka{^) fits 
the experimental curve at R\\ ~ 0. We find that if we 
integrate up to 0.25 eV above the Fermi level we can 
quantitatively explain the data, as seen in Fig. |5| Given 
the limitations of our theory this is essentially equivalent 
to consider ReFpfca(w) for adatom on the Cu(lOO) surface 
for bulk states as a free parameter of the problem. More 
accurate calculations, beyond the ones presented in this 
paper, are required in order to calculate Rc-FpfeaC'-^) from 
a microscopic theory. 
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Figure 14: Differential conductances for the Co on Cu(lOO) for 
different values of . From bottom to the top: -R|| = 0, 2, 4, 6 
A 

In Fig. El we show the amplitude A{R\^) for Co on 

Cu(lOO) for R±^ = 9 A. There is a qualitative agreement 
between theory and the experiment. In Fig. 1141 we show 
the evolution of the tunneling conductance as a function 
of the distance from the adatom. Once again, the line- 
shape is correct for short distances but for long distances 
the lineshape becomes a positive lorentzian, in contrast 
with the experiment. In fact, in Fig. 1151 we show the de- 
pendence of the Fano parameter, q, as a function of the 



lateral distance, R\\ . Notice that overall result is better 
than the case of Co on Cu(lll) fFig. I12|) but the diver- 
gence in the value of q is very clear. 
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Figure 15: Plot of the Fano parameter q{R\\) for Co on 
Cu(lOO) as compared with the values obtained by the ex- 
perimental fili^i. 



VI. CONCLUSIONS 

In this paper we have studied the microscopic theory 
of the single impurity surface Kondo resonance starting 
from an Anderson impurity Hamiltonian. We have calcu- 
lated the hybridization energies and wavefunctions of the 
model adopting simple models for the crystal and surface 
potentials in the nearly-free electron approximation. Our 
model reproduces the main features of the band structure 
of Cu, the energy of the surface states and surface reso- 
nances as measured in photoemission experiments. 

We have shown that the hybridization energies (or ma- 
trix elements), calculated from the solution of the single 
particle Schrodinger equation, are strongly dependent on 
the direction of the electron momenta, the distance of the 
adatom from the surface, and the actual surface where 
the magnetic impurity is located. The hybridization ma- 
trix elements oscillate with distance from the surface due 
to the interference between the adatom localized state 
and the substrate states (surface and bulk) that oscillate 
close to the surface. We have shown that the surface 
Kondo effect occurs preferentially with the surface state 
(when it exists) or with the bulk state at the Fermi en- 
ergy that has the largest component of its momentum 
perpendicular to the surface. 

We have demonstrated that although the surface 
Kondo effect can be obtained from the Anderson im- 
purity Hamiltonian by a Schrieffcr-Wolf transformation 
onto the Kondo Hamiltonian, the hybridization energies 
can be comparable with the adatom atomic energy re- 
quiring a more sophisticated approach that allows for 
fluctuations in the mean occupation of the adatom (in 
other words, the adatom is in what is called a mixed- 
valence state). We have used the mean- field slave-boson 
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formalism that allows for the calculation of the Kendo 
temperature in the limit when the on-site Coulomb en- 
ergy U is much larger than the hybridization energy. Our 
results produce reasonable values of the Kondo temper- 
ature for these systems and the distance of the adatom 
from the surface. 

We have calculated the tunneling conductance of a 
STM tip for different surfaces in the presence of a Co im- 
purity and found that for the surface state in the Cu(lll) 
surface we get quantitative agreement without any fur- 
ther adjustable parameters. However, for the adatom on 
the Cu(lOO) surface where the Kondo effect is dominated 
by the bulk state with momentum along the (100) direc- 
tion, we have to add an extra cut-off energy in order to 
fit the data. More detailed band structure calculations 
are required in order to eliminate the need for this extra 
parameter. We have shown that in both cases we can 
fit very well the data obtained in STM measurements at 
= but there is a discrepancy between the theory and 
experiments when the distance between the tip and the 
adatom becomes large. While in the experimental paper 
where the measurements are reported it is conjectured 
that the Kondo effect in the Cu(lll) and Cu(lOO) sur- 
faces are dominated by the bulk states^!, we show that in 
fact, the Kondo effect is dominated by the surface state 
in the Cu(lll) surface and by the bulk state with mo- 
mentum in the (100) direction in the (100) surface. We 
show that this occurs because there is very little differ- 
ence between surface states and bulk states close to the 
surface in a noble metal and both states have maximum 
penetration on the work function potential. 

In summary, we provide a microscopic interpretation 
for the STM experiments of surface Kondo effects. Our 
results can be used as basis for the interpretation of more 
complicated situations where, for instance, Co adatoms 
interact via direct or indirect exchange on artificial An- 
derson lattices in metallic surfaces. 

Acknowledgments We would like to thank I. Affleck, 
P. W. Anderson, A. Balatsky, P. Coleman, D. Eigler, 
A. Heinrich, P. Johnson, C. Lutz, V. Madhavan, P. Si- 
mon, and A. Tsvelik, for illuminating discussions. We 



also would like to acknowledge support under DARPA 
contract no. DAAD19-01-C-0060. 



Appendix A: EXPLICIT COPPER 
WAVEFUNCTIONS 

For the Cu(lll) surface state, the parallel momentum 
k|| is a good quantum number, and the wavefunction can 
be solved from potential Eq. jnj as: 

^(r) = e*"-""" 

e'^^ cos{kzZ + 6), for z < 0, 

, ylcos Az -I- ;Bsin Az, for < z < zi, 

Cz'e--'t/(l-4^,2;2«:z'), for z > z^. 

(Al) 

where U{a,b; x) is the standard 2nd solution of confluent 



hypergeometric equatioi 



^, z' 



Zim, and kz = tt/2. 



Inside the crystal one has a phase shift S, a crystal decay 
factor /i, and the energy of the surface band bottom E 
related by 



yiC0s25 



^ ■ sm^ 26 + - E, 



2m 



2^2,^2 



2ma2 



-E 



1/2 



\ 2ma^ 



EU2) 



where m = 0.84rne is the bulk-band effective mass. The 
plane wavefunction has a wavevector 



A = i \2m(E + 2Vi 

h 



11/2 



while vacuum decay factor k is defined as 



K^l[2m {Vo - E)] 



1/2 



(A3) 



(A4) 



The energy E is obtained from the continuity conditions 
at z = and z = zi. 
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